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Abstract. - The recently observed multiferroic behavior in the s—1/2 chain cuprate LiVCu04 
(= LiCuV04) with edge-shared Cu04 plaquettes and helical spin ordering does not agree with 
the existing theories such as a spin-current scenario. We argue that the effect can be consistently 
explained, if the exchange-induced electric polarization on the out-of-chain Cu^"'' centers substi- 
tuting for Li-ions in LiVCu04 is taken into account. These substituent centers are proved to be 
an effective probe of the spin incommensurability and the magnetic field effects. 
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The recent observations of multiferroic behaviour in 
edge-shared chain cuprates LiVCu04 [1,2] and LiCu2 02 
[3] challenge the multiferroic community. At first sight, 
these quantum s=l/2 ID helicoidal magnets seem to be 
prototypical examples of ID spiral-magnetic ferroelectrics 
revealing the " ferroelectricity caused by spin-currents" [4] 
with the textbook expression for the uniform polarization 
induced by a spin spiral with the wave vector Q: 



P cx [ea X Q] , 



(1) 



where 63 is a vector orthogonal to the spin spiral plane [5] . 
However, both LiVCu04 and LiCu202 show up a behavior 
which is obviously counterintuitive within the framework 
of spiral-magnetic ferroelectricity [2,3]. Indeed, the spon- 
taneous ferroelectric polarisation in LiVCu04 is observed 
along the a-axis in agreement with the spin-current sce- 
nario. Pa decreases with increasing external field along 
a-axis and vanishes for h > 2 T, where the 6c-plane heli- 
cal structure is realized, however, without any signatures 
of the ferroelectricity along c-axis predicted by the exist- 
ing theories. The saturated value of Pa in LiVCu04 in 
h = reveals a magnitude (Pa « 43/xC/m^) compara- 
ble with that of the multiferroic Ni3V208 where the mag- 
netic ordering drives the electric polarization « lO^/xC/m^ 
(ref. [6]) that represents a typical magnitude of polar- 
ization induced by magnetic reordering in multiferroics. 
However, such an anomalously strong magnetoelectric ef- 
fect seems to be an unexpexted feature for a system with 
Cg-holes and a nearly perfect highly symmetric chain struc- 




Fig. 1: (Color online). Crystal structure of LiVCu04. (In the 
on-line colour version, the ions are highlighted in red, the 
Li"*" ions in blue, the VO4 tetrahedra are highlighted in green, 
the Cu02 chains in light olive). 



ture (see Fig. 1) which both are unfavourable for a strong 
spin-electric coupling. Thus the magnetoelectric effect in 
the title cuprate LiVCu04 raises fundamental questions 
about its microscopic origin. 

Currently two essentially different spin structures of net 
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electric polarization in crystals are considered: 1) with a 
bilinear nonrelativistic symmetric spin coupling [7-9] 



(2) 



or 2) a bilinear relativistic antisymmetric spin coupling 
[4,5,10] 

p„ = ^[n;;,„x[s„xs„]], (3) 

mn 

respectively. The effective dipole moments Tlf^ depend 
on the m, n orbital states and the mn bonding geom- 
etry. First term stems somehow or other from a spin 
isotropic Heisenberg exchange interaction, see, e.g. ref. 
[11], and is not at work for Cu02 chains in LiVCu04 be- 
cause of inversion symmetry consideration. The second, 
or "spin- current" term is at present frequently consid- 
ered to be one of the leading mechanisms of multifcrroic- 
ity [4,5,10,12-14]. However, there are notable exceptions, 
in particular the manganites RMn205, HoMnOs, where a 
ferroelectric polarization can appear without any indica- 
tion of magnetic chiral symmetry breaking [7, 9] , and de- 
lafossite CuFei_2;Ala:02, where the helimagnetic ordering 
generates a spontaneous electric polarization || to the he- 
lical axis [15], in sharp contrast with the prediction of the 
spin current model. At closer examination the spin current 
mechanism comes to a dipole polarization of the three- 
center M1-O-M2 system due to an exchange-relativistic 
antisymmetric Dzyaloshinsky-Moriya (DM) type coupling 
with or without lattice degrees of freedom involved [4, 10]. 
However, the original "spin-current" model by Katsura et 
al. [4] and its later versions [13, 14] seem to be questionable 
as the authors proceed with an unrealistic scenario [16]. 
At variance with such a scenario one can apply a con- 
ventional procedure to derive an effective spin-operator 



P„ = 



na 
12 



Si X S2 



for the electric dipole moment 

in the three-center M1-O-M2 system like a technique sug- 
gested in refs. [17, 18] to derive expressions for the Cu and 
O spin-orbital contributions to the Dzyaloshinsky-Moriya 
coupling in cuprates. However, a straightforward applica- 
tion of this technique to the system of the e^-holes in a 
perfect chain structure of edge-shared CUO4 plaquettes as 
in LiVCu04 shows that the spin current does not produce 
an electric polarization [16]. An alternative mechanism of 
giant magnetoelectricity based on the antisymmetric DM 
type magnetoelastic coupling was proposed recently [10]. 
However, here we meet with a "weak" contributor. Indeed, 
the minimal value of 7 parameter (7 = dD/dH) needed to 
explain the experimental phase transition in multiferroic 
manganites is two orders of magnitude larger than reason- 
able microscopic estimates [10]. Thus we may state that 
the ideal edge-shared Cu04 plaquettes chain structure ap- 
pears to be robust regarding the proper spin-induced elec- 
tric polarization. It means that we should look for the 
origin of the puzzling multiferroicity in LiVCu04 some- 
where within the out-of-chain stuff. 



In the Letter we argue that the ferroelectricity in this 
compomid is actually unrelated to spin currents. It can 
be consistently explained if one takes into account the 
nonrelativistic effects of the spin-dependent electric po- 
larisation [16] due to the out-of-chain Cu^+ centers sub- 
stituting for Li-ions in LiVCu04 present in real samples. 
For an illustration we start with a short description of a 
simple theory of the exchange-induced electric polarisa- 
tion (2) which generalizes the model approach by Tanabe 
et al. [11]. For simplicity, let us consider an one-particle 
(electron/hole) center in a crystallographically ccntrosym- 
metric position of a magnetic crystal. Then all the parti- 
cle states exhibit a definite spatial parity, even (g) or odd 
(u), respectively. Having in mind the 3d centers, we'll as- 
sume an even-parity ground state \g). For simplicity we 
restrict ourselves by only one excited odd-parity state \u). 
An isotropic exchange coupling with the surrounding spins 
can be written as follows: 



(4) 



where /(Rn) is an orbital operator with a ma- 



trix 



The parity-breaking off- 



^ug (Rn ) -^uu (R71 ) 

diagonal part of exchange coupling can lift the center 
of symmetry and mix \g) and \u) states giving rise to a 
nonzero electric dipole polarization of the ground state 



P = 2cgu{g\er\u) = ^ n„(s • S„) , 



where 



Iln — 21 gu (Rn ) 



{g\er\u) 



(5) 



(6) 



and A„g = e„ - gg. 

It should be noted that at variance with the DM type, 
or the spin current scenario the direction of the exchange- 
induced dipole moment for i,i pair does not depend on 
the direction of the spins Si and Sj . In other words, the 
spin-correlation factor (Sj • Sj) modulates a pre-existing 
dipole moment Iljj which direction and value depend on 
the Mej-O-Mcj bond geometry and the orbitals involved 
in the exchange coupling. The net exchange induced po- 
larization of the magnetic crystal depends both on the 
crystal symmetry and the spin structure. 

The orthorhombicaly distorted inverse spinel structure 
of LiVCu04 [19] contains chains of edge-shared LiOe and 
CuOe octahedra running along a and b crystal axes, re- 
spectively [20]. The two neighboring CUO2 chains within 
the same ah plane arc connected by VO4 tctrahcdra that 
alternate up and down along the chain direction, while 
two Cu02 chains within the same unit cell are connected 
by LiOg octahedra (see Fig. 1) First neutron-diffraction 
measurements [22] and a more detailed INS study [23] 
have revealed a long-range incommensurate helical mag- 
netic ordering with a propagation vector q = (0,0.532,0) 
and moments confined to the a6-plane below Tm ~ 2.1 
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K., i.e. the spin spiral runs along the 5-axis: S{y) = 
S'(cos((?), sin(6'), 0) is the classical ground state of an iso- 
lated Heisenberg chain in LiVCu04. Here 9 = qy + a, a 
is a phase shift. The two Cu02 chains within the same 
unit cell are ordered antiferromagnetically: i.e. the phase 
shift Aa = TT. A combined NMR and AFMR study [24] 
also shows that in the low-field range H < Hd « 2.5 T, 
an a6-oriented planar spiral spin structure is realized in 
agreement with refs. 23,24. Based on NMR spectra sim- 
ulations [24] and neutron scattering studies [2], the tran- 
sition at Hci can be attributed to a spin-flop transition, 
where the spin plane of the magnetically ordered struc- 
ture rotates to be perpendicular to the direction of the 
magnetic field applied in the easy plane. 

Noteworthy, LiVCu04 single crystals can reveal a size- 
able non-stoichiometry with Li deficiency and appearance 
of Cu^+ or Cu^+ ions on native Li (Ad) positions [20] due 
to similar ionic radii of Cu and Li. Indeed, a valence bond 
analysis [21] for such a Cu center with unrelaxed nearest 
neighbor O distances of ri—2.l A and r2=2.18 A yields a 
valency i^cu = J2 scu-o ~ 2: 

i^cu = 4exp +2cxp = 1-798, 

where ro{Cu^+)^1.679 A and B = 0.37A. Adopting 
7-1=2.05 A for the Cu-0 distances near the basal plane 
after a small relaxation of the adjacent unshared O ions 
above/below the V ions towards the Cu center one ob- 
tains i/c'„=1.98. These numbers should be compared with 
vcu = 2.16 obtained for a regular Cu chain site. 

Significant, up to 5% concentration of Cu ions on Li- 
sites was assumed by Kegler et al. [25] to explain the ap- 
pearance of an additional weak ^Li NMR signal. The Li 
ions reside in two types of the centrosynimetrical positions 
A,B that differ by a bc-planc mirroring (see Fig. 1). Each 
Li ion has four nearest Cu^+ ions at the same distance. 
Below we argue that the Cu ions substituting for Li posi- 
tions form strongly polarizable entities whose electric po- 
larization due to a parity-breaking exchange interaction 
with the spin spirals in the chains is responsible for the 
multiferroicity in LiVCu04. 

Fig. 2 shows both crystalline surrounding of the Cu^"'' 
substituent center and oxygen-centered hole density dis- 
tribution for active even-parity 629 (oc dxy) and odd- 
parity eu(c) orbitals within the T)ih point group. Let 
us consider the exchange-induced electric polarisation ef- 
fects for these centers starting from a somewhat idealized 
"purely oxygen" model with active 02p 629 and e„((T) or- 
bitals exchange-coupled with two Cu02 chains (see Fig. 2). 
These orbitals are coupled by a large electric dipole matrix 
element: 

{b2g\e:K.\euy{cr)) = (&2g|ey|e„:r(cr)) = --^ei?c«o , (7) 

where RcuO ~ 2 A is the Cu-0 distance. Main exchange 
hole-hole coupling is determined by a "vertical" pdvr trans- 
fer between 02px,y orbitals of central plaquette and Cu 




Fig. 2: (Color online). A slightly idealized view of the Li-0- 
Cu bonds in the unit cell of LiVCu04. An impurity center 
with a Cu-ion substituted for Li is inbetween upper and lower 
Cu02 chains from the same unit cell. Top picture illustrates 
the hole density distributions for active even- and odd-parity 
O 2p orbitals. Main hole p-d transfer channel is shown. 

'idxz,yz (cg-) orbitals of Cu ions in upper and lower chains. 
Diagonal Igg = /&2g!)2g and off-diagonal Igu = higc^ ex- 
change integrals for e.g. Cu^'''-A bond can be written as 
follows: 

Igg{CUn) = tp^^, ( 1 , 

Iguxi^) = Iguyi^) = (8) 

V L + ^ 

where = +74-hB-h2C, A* = -fA-5Barethe 
energies of the two-hole copper singlet and triplet terms 
with 62363 configurations, respectively. A, B, and C are 
Racah parameters, A„g = e^^^ — e^,^^ . The exchange chan- 
nel we consider results in a significant ferromagnetic gg- 
coupling between Cu^+-substituent and adjacent CUO4- 
centers from lower and upper chains. However, the mag- 
nitude of the off-diagonal gw-exchange integrals can suf- 
ficiently exceed that of a conventional diagonal exchange 
integral mainly due to a smaller value of the energy sepa- 
ration (Ag|^* — Aug) as compared with A^^^*. Simple symme- 
try considerations point to the following relations between 
exchange-dipole parameters n„ 

Hi, = H2. = -Ha. = -n4. = ±H ; 
Hiy = -H2g = -Haj, = H4g = H , (9) 
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Fig. 3: (Color online). Illustration of the exchange-induced 
polarisation effect for a Cu 3d 02p h2g hole at a Li-position of 
LiVCu04. The chain spiral Cu^"'' spins are arranged to give 



where upper and lower signs correspond to positions A and 
B, respectively. Illustration of the exchange-induced po- 
larization effect for 02p hole in Li-position of LiVCu04 
is presented in Fig. 3. It should be emphasized that the in- 
duced electric dipole moment is of purely electronic nature 
and lies in the ab plane. To describe different configura- 
tions of the spin neighbourhood for a Cu^+ substituent 
(see Fig. 2), we introduce four basic vectors similarly to 
conventional ferro- and antiferromagnetic vectors as fol- 
lows: 

¥(y) = [Si + S2 + S3 + S4]; G{y) = [Si - S2 + S3 - S4]; 

A(y) = [Si + S2 - S3 - S4]; C{y) = [Si - S2 - S3 + S4] 

(y points to the substituent coordinate along a spin chain) 
with a kinematic constraint: (F • A)=(C • G)=0 valid for 
two identical spirals in the unit cell irrespective of phase 
shift. In terms of new vectors we arrive at simple expres- 
sions for local electric polarization in A and B positions: 



Pa{y) = ±da(s • A) ; Pb{y) = d„(s • C) . 



(10) 



For two ferromagnetically ordered spirals (phase shift 
Aa = 0) A = G = 0, while for two antiferromagneti- 
cally ordered spirals as in LiVCu04 (phase shift Aa — it) 
F = C = 0. It means that the 6-component of the 



electric polarisation Pb{y) for a mutual antiferromagnetic 
(antiphase) spin ordering of the "upper" and the "lower" 
Cu02 chains (Si = — S4,S2 = —S3) vanishes: Pb{y) = 
irrespective of the spin-spiral plane orientation, and the 
substituent electric polarisation may be oriented only 
along crystal a-axis. Moreover, the nonzero net crystal 
polarisation {Pa{y)) one obtains only, if ^ s^, or for 
magnetically nonequivalent A and B substituent positions. 
Spin polarisation of the Cu^^ substituent spin can be eas- 
ily found within the framework of a weak coupling ap- 
proximation, if one takes the most general form of the 
impurity-spiral ground state (gg) exchange interaction 



E 

i=l-i 



Si (i)S, = (s-Ho)^ 



(11) 



where Hq is an effective magnetic field, acting on the 

Cu2+ substituent, laaii) = Iaa,Ixz{i) = ±Ixz\Ixy{'^) = 
~hy{2) = 4^(3) = -hy{A) = ±Ixv;hy{l) = -hy{2) = 

Izy{^) = ~Izy{^) — Izy form a symmetric matrix of the 
exchange integrals, the upper and lower signs correspond 
to positions A and B, respectively. Thus for effective field 
we obtain 

Ho(2/) = IfF -I- TgG (12) 

with 
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Effective molecular field on the Cu^+ substituent spin in 
LiVCu04 is determined solely by anisotropic terms due to 
a complete cancellation of isotropic contributions (F = 0) . 
Thus, the exchange-induced electric polarisation Pa{y) for 
02p b2g hole in Li-position can be written as follows: 

Pa{y) - ±da{s-A) - ± A:^(h-|-Ho)-A = ±^{h-A+ 



'ii{y) 



Hiy) 



'^IxyiGxAy + GyAx) + Izy 

{G,Ay + GyA,))- (13) 

where da = V2ei?c«oA;;g^/bjge„,(<T) (Cwi), h is an exter- 
nal magnetic field. It is worth noting that the effective 
magnetic field in these expressions reflects the symmetry 
of conventional {gg) exchange coupling, while the basic 
spin vectors A and C do that of the parity-breaking {ug) 
exchange coupling. 

We start with an afo-planar spiral as the zero-field 
ground state of the CUO2 chains in LiVCu04 and as- 
sume T=0. Then the nonzero basic spin vectors and the 
resultant magnetic field acting on the spin of the Cu^+ 
substituent at y-position can be written as follows: 

Gx = -45" sin(y ) sin(gy); Gy = AS sin(y ) cos(gy); 
Ax = 4S' cos( y ) cos(gy) ; Ay = 45* cos( y ) sin(gy) . (14) 



Hx = hx + Hox = hx ±ASIxySm{^)cos{qy); 
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Fig. 4: (Color online). The dependence of the mean value of the 
a-axis electric polarisation of the Cu^"*" substituent {Paiy)) (in 
units of da/ion) on the ratio of exchange parameters Ixy/Izy- 
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Fig. 5: (Color online). The field dependence of the mean 
value of the a-axis electric polarisation of the Cu^"'' substituent 
{Pa{y)) (in units of da/ ion) for ab-plane spiral: Ixy/Izy ~ 0.5. 



Hy = hy + Hoy = ky^ I xy Slll ( ) Slll ( ) ; 

H^^h, + Ho, =h,+ 4:SI,y sin( y) cos(gy) , (15) 

We see that both all the components and absolute value of 
effective magnetic field oscillate along the chain direction 
in a complex way thus providing at variance with the DM 
type mechanism the oscillatory behavior of electric polar- 
isation accompanying the spin spiral ordering. Electric 
polarisation Pa (y) oscillates along the chain direction b as 
follows: 

„ ^ ^ , 2daScosif)^^ ^ ^ ^ . , ^ 
Paiy) = ± [hxCOs{qy) + hySm{qy) 

±2I,ySm{^)cos{2qy)]. (16) 

Interestingly, the change of the spiral modulation wave 
length, inversely proportional to the pitch angle (qb), 
should affect the magnitude of the electric polarisation. 
Formally, the spontaneous (h — 0) zero-temperature po- 
larisation Pa{y) does depend on the Ixy/\Izy\ ratio rather 
than on the numerical values of exchange anisotropy pa- 
rameters. It changes sign with that of I^y At the same 
time, (Paiy)) — >0 given I^y — *-0. In Fig. 4 we presented 
the dependence of the mean value of the a-axis electric 
polarisation of the Cu^"*" substituent {Pa{y)) (in units of 
da/ion) given qb/2 = 0.2337r that corresponds to a pitch 
angle ~ 84° [22,23]) on the Ixy/Izy ratio that points to a 
significance of both integrals in the magnetoelectric effect. 

We see that for a reasonable ratio between the exchange 
integrals Ixy and Izy the averaging procedure can signifi- 
cantly reduce the local polarisation from several times to 
more than one order of magnitude. Maximal magnitudes 
of polarisation one obtains for \Ixy/Izy\ ~ 0.4. Fig. 5 
shows the field dependence of the mean value of the a- 
axis electric polarisation of the Cu^"*" substituent {Pa{y)) 
in frames of a rigid a6-plane spiral approximation given 
Ixy/Izy = 0.5. In other words, we assume the external 
magnetic field does not distort the spin spiral, however, 



polarises the substituent spin. It is worth noting that the 
averaging procedure should be performed carefully, as we 
deal with a very slow convergence of the chain summing. 

For the 6c-plane spin spiral ordering the electric polari- 
sation Pa (y) oscillates both for A and B substituent posi- 
tions with a magnitude which can be easily obtained from 
(16), if one makes the interchange:/ia; h,, ±Ixy <-> Izy 
In other words, it means that in a zero magnetic field 
Pailxy/Izy)\ab = Paihy / Ixy)\bc- Howcver, at variance 
with the a6-plane spiral, here we deal with an exact com- 
pensation of A- and B-contributions to the net electric 
polarisation {Pa{y))- {Pa{y)) a = {Pa{y)) b = 0- 

For the ac-plane spiral the molecular field Hg on the 
copper substituent is directed along the b-axis, hence, in 
accordance with (13) and in absence of external magnetic 
field Pa{y) vanishes. The external magnetic field applied 
in ac-plane induces the oscillatory electric polarisation, 
however, the resultant polarisation {Pa{y)) vanishes irre- 
spective of the field direction. 

Thus the spin-spiral ordering in the Cu02 chains of 
LiVCu04 does induce due to the parity-breaking exchange 
a spontaneous electric polarisation on the Cu^+ centers 
substituting for Li-ions resulting in a nonzero net ferroelec- 
tric polarisation P || a-axis only for afc-planar spin spirals, 
and a zero net effect both for ac- and 6c-planar spin spirals. 
Such a behavior differs sharply from the predictions based 
on the spin current scenario (see Exp.(l)), but it is strictly 
confirmed in experiments [1,2] where the ferroelectric po- 
larisation in LiVCu04 has been observed only along a-axis 
irrespective of the direction of the external magnetic field. 
Our model predicts the spontaneous electric polarisation 
Pa is proportional to the m?{T)^ m{T) being the mag- 
netic order parameter, which agrees with the experimental 
findings [2]. The model reproduces nicely a smooth sup- 
pression of the effect with increasing external field h || c. 
Given reasonable estimates for the off-diagonal exchange 
integrals: Ib^ e^{<T) ~ 0.1 eV, and for the g — u energy 
separation: A„g « 4 eV, we arrive at an estimate for the 
maximal value of the electric polarisation per unit cell oc- 
cupied by an impurity P « O.leA, or P « 0.5 /cm? [26]. 
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However, the actual net polarisation is strongly reduced 
compared with this huge magnitude due to the small im- 
purity concentration c <C 1 and the oscillatory character 
of Pa{y)- Anyhow, even for reasonable small substitutions 
c(Cuii)ss 0.03 — 0.05 [20,25] and after averaging over im- 
purity positions the net impurity polarisation in LiVCu04 
can still be comparable with that observed in the gigantic 
multiferroic Ni3V208 [6]. 

Above we assumed a weak coupling between the Cu^"*" 
substituent and the adjacent spin spirals treating the sub- 
stitucnt as a paramagnetic impurity. However, the Cu^"*" 
substituent itself can locally distort the spin spirals. Irre- 
spective of the sign of the exchange parameters the Cu^+ 
substituent - spin spiral exchange interaction favours lo- 
cally a ferromagnetic inter-chain coupling with a phase 
shift Aa = 0, thus leading to a inter-chain coupling com- 
peting with that of the ideal chain system. The role played 
by this effect in magnetoelectric coupling deserves a spe- 
cial investigation. Also the effect of an external magnetic 
field on the ferroclcctricity in spiral magnets seems to be 
more involved than it appears in our simplified model the- 
ory, and needs both in a more detailed and sophisticated 
description of interacting impurity centers, spirals, and a 
field induced changes of the spirals themselves. In partic- 
ular, the pecularities near 8T observed in the magnetiza- 
tion measurements [27,28] may evidence a breaking of the 
weak Cu^+ substituent - spin spiral coupling approxima- 
tion thus making a breakdown of the ferroelectricity ob- 
served near ST [1,2] quite natural. Finally, also the effect 
of finite temperatures should be considered. As regards 
the quantitative estimates of the net electric polarisation 
in LiVCu04 it is worth noting that, strictly speaking, the 
parity-breaking exchange results in electric polarisation of 
both exchange-coupled centers with a trend toward a can- 
cellation, if the inversion center takes place inbetwcen. In- 
deed, the Cu^+ substituent will polarise the four adjacent 
Cu04 centers which net polarisation, however, fails to can- 
cel Pa (y) due to a nonequivalence of the exchange-coupled 
Cu04 centers. 

Anyhow, we have shown that in contrast to a spin- 
current scenario our simplified microscopic model already 
catches qualitatively correct the main features of the mul- 
tiferroicity observed in LiVCu04, namely the a-axis direc- 
tion of the net electric polarisation for the a5-planc spin 
spiral ordering and its vanishing for the be- and ac-planar 
spirals. Along with a reasonable estimation of the po- 
larization, the model is believed to explain the subtleties 
of the field dependence, in particular, the suppression of 
the polarisation with increasing external magnetic field || 
c-axis. 

In conclusion, we have discussed recent observations 
of the multiferroic behavior in the ID chain cuprate 
LiVCu04 with the edge-shared CUO4 plaquettcs and in- 
commensurate spin ordering and argued that its peculiar- 
ities can be consistently explained if one takes into ac- 
count the exchange-induced electric polarisation on the 
Cu^+ centers substituting for Li- ions in LiVCu04. Thus 



we deal with an unusual system where in sharp contrast 
with manganites the multifcrroicity happens in two dif- 
ferent subsystems, that is magnetic spirals induce elec- 
tric polarisation at Cu^+ ions substituting for native Li+ 
ions. These substituent centers arc proved to be an effec- 
tive probe of the spin incommensurability and magnetic 
field effects. The experimental investigation of possible 
relations between physical properties and the extent of 
nonstoichiometry as well as the theoretical consideration 
of analogous scenarios in other multiferroics would be of 
considerable scientific and potential practical interest. 

Note added. After the submission of the manuscript, we 
became aware of an independent experimental study of the 
multifcrroicity in LiVCu04 by Schrettle et al. [29] . In con- 
trast with Yasui et al. [2] , the authors have found a weak 
c-axis electric polarization also for a 6c-plane spin spiral. 
Generally speaking, this discrepancy does not necessarily 
contradict the nonstoichiometric nature of the multifer- 
roic behaviour which makes it strongly sample dependent. 
In the general framework of an impurity model this might 
point to a manifestation of weak charge transfer effects be- 
tween the Cu^"*" substituent and the spin spirals. Anyhow, 
further theoretical and experimental studies are needed to 
uncover all the subtle details of the unusual and challeng- 
ing multiferroicity in LiVCu04. 
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